Improved detection of small atom numbers through image processing 
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We demonstrate improved detection of small trapped atomic ensembles through advanced post- 
processing and optimal analysis of absorption images. A fringe removal algorithm reduces imaging 
noise to the fundamental photon-shot-noise level and proves beneficial even in the absence of fringes. 
A maximum-likelihood estimator is then derived for optimal atom-number estimation and is applied 
to real experimental data to measure the population differences and intrinsic atom shot-noise be- 
tween spatially separated ensembles each comprising between 10 and 2000 atoms. The combined 
techniques improve our signal-to-noise by a factor of 3, to a minimum resolvable population differ- 
ence of 17 atoms, close to our ultimate detection limit. 



Trapped ultracold atoms and quantum degenerate 
gases are novel systems for the study of many-body 
quantum physics 1] and are key to new technologies 
such as trapped atom interferometers [_2!| and atomic 
clocks 3]. This is exemplified by recent experiments on 
number/spin-squeezing and entanglement between small 
atomic ensembles [J| , which could serve as a resource for 
quantum metrology and quantum information science. 
Essential to such applications is the ability to precisely 
measure the populations of a pair of atomic ensembles 
in, for example, a double- well potential for Josephson 
physics or atom interferometry [2|, mW, or the spin states 
of an atomic clock. Upon readout the interferometric 
phase can be mapped to a population difference and 
the two ensembles are imaged to different regions of a 
CCD camera. The relative population is robust against 
common-mode technical fluctuations such as probe noise 
or trap loading efRciency. However, imaging noise and in- 
trinsic atom number fluctuations (atom shot noise) typ- 
ically limit measurement precision. 

Here we demonstrate improved detection of trapped 
ensembles of ultracold atoms through advanced post- 
processing and optimal analysis of laser-illuminated ab- 
sorption images. This allows us to better measure the 
intrinsic atom number fluctuations in a magnetic lattice 
potential. First we apply a fringe-removal algorithm to 
reduce residual imaging noise to the fundamental photon 
shot-noise level. We then establish the ultimate limit for 
measuring the relative populations based on the Cramer- 
Rao bound (CRB) and derive maximum likelihood esti- 
mators used to attain this limit [S] . These optimal anal- 
ysis techniques provide the basis to improve the readout 
of trapped atom interferometers to the quantum limit or 
to better resolve number squeezing and entanglement. 

In our experiment we prepare up to a total of A'^ = 
4 X 10^ ^^Rb atoms in a multi-well trap formed by a cur- 
rent through a Z-shaped wire and the permanent field 
of a magnetic lattice atom chip (9J. Recently we used 
this setup to study sub-Poissonian atom number fluctua- 
tions resulting from three-body loss in an array of tightly 
confining microtraps [lO|. In this work we position the 



12 



10 



s 

O 

e 



FIG. 1. Typical absorption image of the atomic distribution 
(90 X 71 pixels) for a hold time of 1.6 s. Four independent 
atomic ensembles are visible (horizontally distributed) with 
two mirror images (vertically separated). The two center- 
most clouds each contain « 400 atoms. Dashed boxes indicate 
the background pixels for the fringe removal algorithm and 
the signal regions are indicated by solid lines. 



trap at the edge of the magnetic lattice, 45 /im below 
the chip surface, to create four potential wells (two deep 
central wells and two shallow outer wells). The atoms 
are evaporatively cooled in this trap to a temperature of 
~ 1 /zK, thereby creating independent atomic ensembles, 
each occupying an area of ~ 10 x 50 /xm^ in our images 
(fig. [1} . The mean number of atoms is varied by holding 
the atoms in the trap for 0-4 s at a fixed temperature 
and partially by reducing the amount of Rb dispensed 
during loading. 

Absorption imaging is the standard method for de- 
tecting trapped neutral atoms. Here atoms are briefly 
exposed to a nearly homogeneous probe laser, typically 
tuned to resonance with an atomic transition and the re- 
sulting absorption signal A is imaged onto a CCD cam- 
era. Subsequently the atoms are ejected from the trap 
and a reference image R is recorded to normalize inten- 
sity variations of the probe. A dark image may also be 
recorded without the probe to subtract any stray light or 
CCD dark counts. The two-dimensional atomic de nsity 
is calculated as n = —a^og{A/ R) + s{A/R — 1)]/(To [ll 



Here ctq = 3A^/27r is the absorption cross-section, a 
is a dimensionless parameter which depends for exam- 
ple on the probe polarization and s — Iji/al'^ is the 
(spatially dependent) saturation parameter, with /^ the 
probe intensity at the position of the atoms and 1° = 
1.67 mW/cm^ the saturation intensity for our transition. 

The experiments employ a back-illuminated deep- 
depletion CCD camera with a quantum efficiency of 
Qe ~ 0.9 at A = 780 nm, a measured gain oi g — 
0.87 ± 0.05 counts/photon and a readout noise level of 
(Trd ~ 13 counts. The resolution of our optical sys- 
tem is 9.6 /im (Rayleigh criterion) with a pixel area of 
A = (3.3 ± 0.1 /zm)^ in the object plane. The probe is 
slightly inclined with respect to the gold-coated chip sur- 
face to create two mirror images of the atoms (fig. [1]) [12| . 
Correlating these images provides additional means to 
distinguish atom shot noise and imaging noise. 

We have experimentally optimized our imaging param- 
eters to provide the highest signal-to-noise ratios (SNR) 
for in-trap imaging. Long exposure times r tend to lower 
the effect of photon shot noise, however blurring/heating 
due to photon recoil during the imaging pulse tends to 
increase the effective area a over which the atoms are dis- 
tributed; therefore, one finds an optimum intensity and 
exposure time. We measure the SNR in a series of ab- 
sorption images for constant atom number as a function 
of both T and s and find a maximum for r = 50 /iS and 
s ~ 0.54. Finally, we calibrate a by comparing the in- 
tegrated absorption signal with a second set of images 
taken after free expansion using a weak probe (s « 0.02), 
yielding a — 3.0 ± 0.2. Detection linearity with respect 
to atom number variations is confirmed by measuring 
a fixed ratio between the populations of the center and 
outer wells. Our typical atom numbers correspond to 
weak absorption signals {A/R > 0.8) for which it can be 
difficult to extract a signal buried in noise. 

In practice additional imaging noise originates from 
fringes due to diffraction and interference of the probe 
beam by optical elements and the atom chip surface. 
Small vibrations between the absorption and reference 
images result in imperfect normalization giving rise to 
fluctuating fringe patterns. This noise can be greatly 
reduced through the application of a fringe-removal al- 
gorithm, while making no assumptions about the atomic 
distribution. The algorithm works by composing for each 
absorption image a matching optimal reference image Q 
constructed as a linear combination of many reference im- 
ages Rk within a set, Q = J^k ^kRk l^uM,- The method 
is closely related to that applied to facial recognition 15 1 
and recently in astronomical image analysis for detecting 
extrasolar planets [l6|. 

To obtain the coefficients Ck we minimize the least 
squares difference between the absorption and reference 
images ^^rrixiAx — QxY, where x indexes each pixel, 
within a specified background region (m^ = 1) excluding 
the signal region {rrix = 0). Setting partial derivatives 



10" 



O" 



10^., 



10' 



10" 




■ Raw images ''"^^'^v 

• Fringe removal 

Predicted noise 

Ultimate limit 



100 



300 



1000 

(R) 



3000 



10000 



FIG. 2. Imaging noise for various probe intensities with 
and without fringe removal, var(^/(3) (circles) and var(A/i?) 
(squares) respectively. The predicted photon shot-noise plus 
camera readout noise is shown as a solid line. The dashed 
line shows the ultimate photon shot-noise limit. 



with respect to c, to zero, we obtain a linear system of 



: X^a; "fTixRxjAx with the square 
Rx,jRx,k, which can be readily 



equations, J^k'^kBj^k = 
matrix Bj^k = Yjx"^^ 
solved for Ck- A typical data set consists of hundreds 
of absorption images, therefore we decompose B once 
using LU or singular value decomposition then substi- 
tute to obtain Ck for each absorption image 17[. The 



algorithm is sufficiently fast to decompose B and pro- 
cess new absorption images in < 1 s for live processing 
between experimental cycles. 

Figure [2] shows the imaging noise with and without 
fringe removal, var(A/(5) and var(A/i?) respectively, cal- 
culated for a signal-free region (separate from the fringe- 
removal background region) for various probe intensities. 
Without fringe removal the measured noise is in good 
agreement with the expected photon shot noise and read- 
out noise given our CCD parameters. Application of the 
fringe removal algorithm with a basis of ~ 250 refer- 
ence images reduces the measured variances by a factor 
of 1.9 ± 0.3 over the full range of intensities. Remark- 
ably, even in the absence of fringes, the algorithm reduces 
the photon shot-noise contribution originating from R. 
This is possible since the optimal reference image is the 
(weighted) average over many reference images, allowing 
an additional decrease of uncorrelated noise by up to a 
factor of 2. For our chosen imaging parameters the re- 
maining noise from pixel to pixel is cr„ — 1.3 ±0.3 atoms, 
close to the ultimate limit of 1.1 atoms due to photon 
shot noise in A. There is also a small residual corre- 
lated noise component which fluctuates on a length scale 
comparable to our cloud size with an rnis amplitude of 
0.06 ± 0.01 atoms. 

We wish to measure the total and difference popula- 
tions N^ ~ Np ± Nq and the relative population differ- 
ence N~ /N^. To establish the limit for extracting these 
populations from our images we derive the Cramer-Rao 



bound (CRB). The CRB is a powerful tool which gives a 
lower bound for the variance of any parameter estimate, 
independent of the exact procedure used to extract the 
information g]. We describe our image data as the ab- 
sorption signal from a pair of atomic ensembles corrupted 
by uncorrelated additive gaussian noise di^x with variance 
<^n'- ''^i,x = Ni^pPx + Ni^qQx + di^x, wherc i labels a partic- 
ular realization. The normalized spatial mode functions 
Px and qx {^xPx = Hx 1^ = 1) are determined by the 
cloud shape and optical resolution of our imaging system. 
The log-likelihood function is 
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The CRB is then the inverse of the Fisher information 
matrix, calculated from the second derivatives of I with 
respect to N^. We obtain for the CRB 
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with the parameters u* = ^^(pa:±9x) andw — X]x(Px 
ql). Here var(Af+) > Cn, var(iV-) > C22, and for the 
relative atom number. 
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Taking typical numbers from our experiment (v ~ and 
l/u"^ « 30) we find the single-shot CRB to be Cn = 
C22 — 200 atoms^, or a minimum resolvable population 
difference of 14 atoms. The CRB is also easily applied to 
estimating the number of atoms within a single ensemble 
for which our detection limit is 10 atoms/shot. 

The CRB can be attained by a maximum likelihood 
estimator (MLE), found by maximizing the log-likelihood 
function and solving for N , 
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(4) 
Equation @ can be interpreted as a sum over the im- 
aged density distribution weighted by Px^qxi where the 
second term accounts for overlapping and uneven mode 
functions. This gives a direct measure of the popula- 
tions while minimizing the influence of noise. The mode 
functions px and qx can be obtained from the data in a 
model independent way by averaging over many images 
to suppress noise, lipx and qx are not spatially separated 
one could record a set of images where each ensemble is 
individually populated or apply independent component 
analysis to isolate signal components [18|. 

We compare the expected performance of the MLE to 
naive estimates for N obtained by separately integrat- 
ing over rectangular sub-images containing the left and 
right ensembles. For this case, the expected variance is 
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FIG. 3. Relative variance as a function of total atom num- 
ber measured by maximum-likelihood estimation (squares) 
or straight integration (circles). The dash-dotted line shows 
the expected fluctuations combining atom shot-noise (dashed 
line) and the minimum detection noise contribution given by 
the CRB (dotted line). Solid lines indicate fits for the detec- 
tion noise contributions. The inset shows the covariance be- 
tween A''" /N'^ in the upper and lower mirror image (squares) 
along with the predicted atom shot noise (dashed line). 



var(iV^) > aa^, with a the total number of pixels in 
the integration regions. For regions chosen to include 
> 95% of the atomic signal (a/2 = 10 x 22 pixels, fig. [T]) 
we expect a detection noise contribution to var(A^^) of 
730 atoms^, a factor 3.7 larger than the CRB. 

We have performed measurements of the relative atom 
number N~ /N~^ for varying total atom number to resolve 
intrinsic atom number fluctuations and demonstrate the 
improvements gained using maximum-likelihood estima- 
tion. Our data set spans 40 hold times between 0-4 s, 
corresponding to varying {N~^) from 4 x 10'^ down to ^ 20 
atoms. For each hold time we repeat the experiment 40 
times. Sub-pixel image registration is applied to align 
each image to the average to eliminate small fluctuations 
in the atom cloud positions. The mode functions are 
computed by averaging all the data and segmenting the 
result. We then extract from the images the atom num- 
ber populations using eq. ^ and compute the variance 
vai{N~ /N~^). The expected atom shot-noise contribu- 
tion to the variance, accounting for a mean population 
imbalance, is 1/(A^+) - {N-)^/{N+y. 

Figure [3] shows the measured variance of the rela- 
tive population for both the MLE and simple integra- 
tion as a function of (iV+). We restrict our analysis 
to the lower mirror image specified in fig. [1] For large 
{N'^) both estimation methods yield results which are 
atom shot-noise limited. For the integration method 
the measured variance is limited by detection noise for 
(N^) below ~ 1000 atoms. Detection is significantly 
improved using the MLE enabling atom-shot-noise lim- 
ited detection down to below 300 atoms. To the data 
we fit the model var(iV-/7V+) = C22/{N+y + 1/{N+) 



(for our data {N ) sa 0), to measure the detection 
noise contribution. For straight integration we find 
(5*22 = (1-3 ± 0.2) X 10^ atoms^, whereas for the MLE 
C22 = 270 ± 40 atoms^, corresponding to a minimum 
resolvable difference of 17 atoms/shot. The measured 
noise is slightly larger than the minimum value given by 
the CRB, but is consistent with the added effect of the 
residual correlated noise. 

We can further distinguish atom number fluctuations 
from detection noise in our data by comparing the mea- 
sured population differences in the two mirror images. 
Here we expect the detection noise to be uncorrelated, 
and true atom number fluctuations to be correlated be- 
tween the two images. To verify this, we measure N~^ 
and N^ for both mirror images separately and calculate 
the covariance of N~ /N~^ between the upper and lower 
images for each hold time. The result is shown in the 
inset of fig. [U along with the expected atom shot-noise, 
l/(iV+). Our data is consistent with atom shot noise over 
the full range of atom numbers, demonstrating a robust 
method to observe intrinsic atom number fluctuations. 

To summarize, we have demonstrated improved detec- 
tion of small atomic ensembles, to close to the ultimate 
photon shot- noise limit. A fringe removal algorithm and 
maximum likelihood estimation are applied to absorp- 
tion images to reach a measured detection sensitivity 
of 17 atoms/shot for population differences. The mea- 
sured variance of the relative populations is in excellent 
agreement with the lower limit given by atom shot noise 
and the Cramer-Rao bound for our imaging parameters. 
Combining fringe removal and maximum likelihood esti- 
mation we have improved our signal-to-noise ratio by a 
factor of 3 allowing for atom shot-noise limited detection 
of ensembles comprising as few as 270 atoms, a factor 9 
lower than without these methods. Averaging measure- 
ments from both mirror images would offer a further y^ 
improvement in SNR. 

The sensitivity could be improved to the single-atom 
regime by increasing the imaging resolution and better 
localizing the atoms. Optimal imaging parameters are 
found by considering a model for blurring due to pho- 
ton recoil after release from the trap, where the cloud 
size a increases proportional to st^/{1 + s). Integrated 
photon shot-noise scales with ■y/a(l + s)^ /st, yielding an 
optimal saturation parameter of Sopt — 2/3. For real- 
istic imaging parameters (optical resolution of 1.2 /im, 
Qe = 0.9 and a = 1) we find an optimal exposure time of 
13 /iS and a detection limit of 0.5 atoms/shot for a single 
ensemble. Further improvements would be possible using 
squeezed light to image ultracold atoms [l9|. 

The atom clouds we prepare are considerably elon- 
gated, highlighting the potential application of our anal- 
ysis to the study of ID quantum gases. Applying maxi- 
mum likelihood estimation column-by-column provides a 
robust measure of linear density. For our experimental 
conditions we infer a sensitivity of 0.8 atoms//im. This 



is below the typical linear density required to reach the 
crossover from the weakly interac ting to the strongly in- 
teracting regime on an atom chip [20| , which could be di- 
rectly imaged in a single realization. By applying a weak 
optical lattice along the length of the cloud, it would be 
possible to directly observe the ID Mott insulator phase 
transition [ij, l2ll | , with a sensitivity of a single atom per 
lattice site. 
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